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When  laser  light  propagates  through  a  resonant  medium,  the  transmitted  beam  can  exhibit  excess  intensity 
noise  [amplitude  modulation  (AM)].  In  a  semiclassical  description  of  the  phenomenon,  laser  phase  noise  (PM) 
induces  fluctuations  in  the  medium’s  electric  susceptibility,  which  in  turn  cause  fluctuations  in  the  transmitted 
intensity.  The  process  provides  an  efficient  means  for  PM-to-AM  conversion,  and  intuition  suggests  that  large 
linewidth  lasers  should  exhibit  much  greater  PM-to-AM  conversion  than  narrow  linewidth  lasers.  Here  we 
measure  the  relative  intensity  noise  (RJN)  for  two  diode  lasers  whose  linewidths  AvL  differ  by  more  than  102, 
after  the  lasers  have  propagated  through  a  resonant  rubidium  vapor.  Though  the  RIN  of  the  narrow  linewidth 
laser  is  only  reduced  by  a  factor  of  about  6  compared  to  the  broad  linewidth  laser,  our  results  are  nonetheless 
consistent  with  numerical  simulations  of  the  PM-to-AM  conversion  process.  In  particular,  both  computation 
and  analytical  theory  indicate  that  RIN  is  a  nonlinear  function  of  A  vL .  For  single-mode  laser  linewidths  less 
than  the  atomic  dephasing  rate,  RIN  increases  like  f/S  vL,  while  for  linewidths  greater  than  the  atomic  dephas¬ 
ing  rate  RIN  is  a  decreasing  function  of  Ai>L.  [SI 050-2947(99)0090 1-4] 


PACS  number(s):  42.50.Gy,  42.62.Fi,  42.25.Bs 

L  INTRODUCTION 

It  is  now  fairly  well  established  that  when  laser  light 
propagates  through  a  resonant  medium,  the  transmitted  beam 
will  exhibit  some  level  of  excess  intensity  noise.  Though  the 
phenomenon  may  be  modeled  in  different  ways,  depending 
on  the  specifics  of  the  experimental  situation,  in  all  cases  the 
excess  intensity  noise  results  from  a  mapping  of  atomic  fluc¬ 
tuations  onto  the  resonant  optical  field.  In  strong  fields, 
where  the  transition  may  be  saturated,  this  phenomenon  is 
described  using  a  quantum  electrodynamics  formalism,  and 
is  explained  as  “vacuum  side-mode  amplification  due  to  for¬ 
ward  four- wave  mixing”  [1].  In  this  regime,  the  phenom¬ 
enon  can  be  quite  significant  even  for  perfectly  monochro¬ 
matic  fields,  as  a  consequence  of  the  medium’s  natural 
relaxation  processes,  and  it  plays  a  role  in  certain  amplitude¬ 
squeezing  experiments  [2].  In  weak  fields  the  phenomenon  is 
no  less  important,  though  it  is  primarily  associated  with  the 
laser’s  intrinsic  phase  fluctuations,  and  can  be  explained 
semiclassically:  laser  phase  noise  (PM)  induces  fluctuations 
in  the  medium’s  electric  susceptibility,  which  in  turn  cause 
fluctuations  in  the  laser’s  transmitted  intensity  [amplitude 
modulation  (AM)]  [3].  At  high  Fourier  frequencies,  the 
weak-field  manifestation  of  the  phenomenon  provides  spec¬ 
troscopic  information  on  the  resonant  medium’s  energy  level 
structure,  as  first  reported  by  Yabuzaki  et  al  [4].  In  optically 
thick  media,  the  low  Fourier  frequency  components  of  the 
excess  intensity  noise  can  be  orders  of  magnitude  larger  than 
the  laser’s  intrinsic  relative  intensity  noise,  and  the  PM- 
to-AM  conversion  phenomenon  becomes  an  important 
physical  process  limiting  the  performance  of  gas-cell  atomic 
clocks  [5,6]. 

The  present  studies  consider  the  role  of  single-mode  laser 
linewidth  in  the  PM-to-AM  conversion  process  for  weak 
fields.  Intuitively,  one  expects  the  magnitude  of  PM-to-AM 
conversion  to  scale  with  linewidth,  since  a  single-mode  la¬ 
ser’s  linewidth  is  principally  related  to  the  laser’s  phase 
noise.  In  particular,  as  suggested  by  the  high  Fourier  fre¬ 


quency  experiments  of  McLean  et  al  [7]  and  McIntyre  et  al 
[8],  if  the  linewidth  of  a  diode  laser  is  narrowed  by  several 
orders  of  magnitude,  one  anticipates  orders  of  magnitude  re¬ 
duction  in  processes  associated  with  PM-to-AM  conversion. 
The  experimental  and  theoretical  results  described  below, 
however,  do  not  completely  bear  out  this  expectation.  In  Sec. 
II  we  describe  an  experiment  that  compares  the  relative  in¬ 
tensity  noise  (RIN)  of  two  diode  lasers  whose  linewidths 
differ  by  two  orders  of  magnitude.  The  measurements  were 
performed  for  several  different  values  of  the  vapor’s  optical 
depth.  On  average,  the  RIN  of  the  narrowed  diode  laser  was 
at  best  only  a  factor  ~6  smaller  than  the  RIN  of  the  much 
broader  diode  laser.  As  discussed  in  Sec.  IV,  these  observa¬ 
tions  are  consistent  with  numerical  simulations  of  the  PM- 
to-AM  conversion  process,  which  additionally  show  that  for 
large  linewidth  fields  RIN  is  a  decreasing  function  of  laser 
linewidth.  The  computational  results  are  elucidated  in  Sec. 
V,  where  it  is  shown  that  the  variance  of  atomic  coherence 
fluctuations,  which  drive  the  PM-to-AM  conversion  process, 
has  a  maximum  when  the  laser  linewidth  equals  the  optical 
transition’s  homogeneous  linewidth. 

H.  EXPERIMENT 

Figure  1  shows  a  block  diagram  of  the  experimental  ar¬ 
rangement.  The  beams  from  a  Mitsubishi  TJS  [9]  AlGaAs 
diode  laser  (ML4102)  and  an  EOSI  [10]  external  cavity  (EC) 
AlGaAs  diode  laser  were  superimposed,  expanded,  and  al¬ 
lowed  to  pass  through  a  natural  Rb  vapor  contained  within  a 
Pyrex  resonance  cell  (72%  Rb85  and  28%  Rb87)  [11].  Only 
one  beam  was  allowed  to  pass  through  the  vapor  at  a  time, 
and  both  lasers  were  tuned  to  the  Rb85  D]  transition  at  794.7 
nm  [i.e.,  5  2Sm(F=  2)— >5  2P\n(F=  2,3)].  As  a  result  of 
Doppler  broadening,  the  excited-state  hyperfine  structure  of 
Rb85  was  not  resolved.  The  absorption  spectrum  of  the  natu¬ 
ral  vapor,  as  measured  by  the  TJS  diode  laser,  is  shown  in 
Fig.  2(a).  The  absorption  feature  of  interest  for  the  present 
work  is  labeled  as  D  in  the  figure,  and  is  well  resolved  from 
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FIG.  1.  Experimental  arrangement  as  described  in  the  text. 

the  other  absorption  features.  An  optical  isolator,  providing 
at  least  36-dB  isolation,  prevented  optical  feedback  from  per¬ 
turbing  the  TJS  laser’s  characteristics,  while  this  function 
was  essentially  performed  by  a  2.0  neutral  density  filter  that 
was  tilted  and  placed  in  front  of  the  external  cavity  diode 
laser.  Prior  to  entering  the  resonance  cell,  a  portion  of  laser 
light  was  detected  with  a  moderately  high  finesse  confocal 
Fabry-Perot  interferometer  (frequency  resolution,  Az>res, 
—  13  MHz).  The  interferometer  allowed  us  to  measure  the 
laser  detuning  from  resonance,  and  also  ensured  that  optical 
feedback  was  not  artificially  broadening  either  laser’s  spec¬ 
trum  during  the  course  of  the  measurements. 

The  TJS  laser  had  a  measured  linewidth  [full  width  at  half 


FIG.  2.  (a)  Absorption  spectrum  of  the  natural  Rb 

vapor  (i.e.,  72%  Rb85  and  28%  Rb87)  as  measured  with  the 
56-MHz  linewidth  diode  laser  at  794.7  nm  and  at  a  temper¬ 
ature  of  32  °C:  (A)  5  2P1/2(F=  l)->5  2S1/2(F=  2)  transition  for 
Rb87;  (B)  5  2P1/2(F=  2)-^5  2S1/2(F=  2)  transition  for  Rb87; 
(O  5  2P1/2(F=2,3)->5  2S1/2(F=3)  transition  for  Rb85;  (D) 
5  2F1/2(F= 2,3)— >5  251/2(F=  2)  transition  for  Rb85;  (F) 

5  2Fi/2(F=1)->5  2S1/2(F=1)  transition  for  Rb87;  (F)  5  2P m(F 
=  2)— >5  251/2(F=1)  transition  for  Rb87.  The  ground-state  hyper- 
fine  splitting  of  Rb85  is  3036  MHz,  while  that  of  Rb87  is  6835  MHz. 
The  5  2P\a  excited-state  hyperfine  splitting  of  Rb85  is  362  MHz, 
and  unresolved  within  the  Doppler-broadened  absorption  line,  while 
that  of  Rb87  is  812  MHz.  (b)  Relevant  energy-level  diagram  for  the 
D  transition  of  (a). 


maximum  (FWHM)],  A  vL ,  of  56  MHz.  Using  a  very  high 
finesse  confocal  interferometer  (A^res,~0.8  MHz)  we  mea¬ 
sured  the  external  cavity  laser  linewidth,  and  found  it  to  be 
limited  by  the  interferometer  resolution.  According  to  the 
external  cavity  laser’s  specifications,  the  linewidth  should  be 
100  kHz  or  less.  Thus  the  two  lasers  had  linewidths  that 
differed  by  at  least  of  factor  of  70,  and  quite  likely  by  nearly 
three  orders  of  magnitude. 

The  laser  power  for  both  beams  was  15  /zW,  measured 
prior  to  the  resonance  cell,  and  both  beams  were  expanded 
and  then  apertured  to  a  diameter  of  0.4  cm.  The  low  light 
intensity  in  the  resonance  cell  ensured  that  neither  laser 
could  cause  any  saturation  of  the  optical  transition.  Addition¬ 
ally,  since  there  was  no  buffer  gas  in  the  resonance  cell,  there 
were  no  confounding  effects  of  optical  pumping  on  PM- 
to-AM  conversion  [12,13].  The  resonance  cell  had  a  length  L 
of  15  cm  and  a  diameter  of  2  cm,  and  was  heated  by  braided 
wire.  For  our  experiments,  the  Doppler-broadened  linewidth 
was  — 510  MHz,  and  the  natural  linewidth  of  the  D\  transi¬ 
tion  is  5.67  MHz  [14]. 

In  order  to  be  sure  that  the  laser  beam’s  passage  through 
the  resonance  cell  had  no  appreciable  effect  on  the  laser 
spectrum,  we  replaced  the  beam  chopper  and  photodiode 
shown  in  Fig.  1  with  one  of  our  two  interferometers.  We  then 
measured  the  linewidth  of  the  TJS  laser  (using  the  interfer¬ 
ometer  with  A  vTQS—  13  MHz)  and  the  EC  laser  (using  the 
interferometer  with  Avres~0.8  MHz)  after  passage  through 
the  vapor.  For  a  vapor  temperature  of  38  °C,  and  for  the 
lasers  tuned  both  off  resonance  and  on  resonance,  the  laser 
linewidths  after  passage  were  indistinguishable  from  the 
linewidths  measured  prior  to  passage  through  the  resonance 
ceH. 

The  vapor  cell  temperature  was  varied  from  22  to  46  °C, 
and  for  these  temperatures  the  transmission  of  the  external 
cavity  laser  was  used  to  estimate  the  vapor’s  optical  depth 
Td :  rj 1  =  L  ~ 1  ln[/ojft7 7on] ,  where  subscripts  refer  to  the  laser 
tuned  either  off  or  on  resonance.  For  this  temperature  range, 
rd!L  varied  from  about  3  to  0.3,  indicating  a  transition  from 
an  optically  thin  vapor  to  an  optically  thick  vapor. 

The  average  laser  intensity  (7)mea s  was  obtained  by  chop¬ 
ping  the  light  beam  and  measuring  the  photodiode  signal 
with  a  lock-in  amplifier.  Then  the  chopper  was  stopped,  and 
the  intensity  noise  SI  in  a  1-Hz  bandwidth  at  1  kHz  was 
measured  with  a  spectrum  analyzer  that  had  a  25-kHz  band¬ 
width.  The  laser  RIN  was  defined  as  the  ratio  of  SI  to  </} 
(where  (I)  corresponds  to  the  measured  average  laser  inten¬ 
sity  corrected  for  amplified  spontaneous  emission  as  dis¬ 
cussed  subsequently),  and  in  all  measurements  SI  was  mea¬ 
surably  larger  than  the  photodetector’s  dark  noise.  (Though 
this  definition  of  the  RIN  is  essentially  the  square  root  of  that 
typically  found  in  the  diode  laser  literature,  in  our  opinion  it 
is  more  amenable  to  spectroscopic  applications.) 

In  order  to  obtain  an  appropriate  value  of  the  RIN  for 
each  of  the  lasers,  it  was  necessary  to  correct  the  TJS  laser’s 
transmitted  intensity  for  amplified  spontaneous  emission 
(ASE).  (This  correction  was  also  performed  for  the  EC  laser, 
but  it  was  of  less  significance.)  As  the  gain  curve  of  a  diode 
laser  is  fairly  broad,  diode  laser  Fabry-Perot  modes  near  the 
main  lasing  mode  will  be  amplified  when  the  laser  is  oper¬ 
ating  above  threshold  [15].  For  the  TJS  laser,  these  modes 
are  spaced  by  about  0.3  nm,  and  so  will  not  be  absorbed  by 
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TABLE  I.  Ratio  of  ASE  to  the  transmitted  intensity  for  the 
lasers  tuned  on  resonance.  [This  relatively  large  value  is  simply  due 
to  strong  absorption  of  the  EC  laser  at  46  °C  (i.e.,  UtTd=  3.8).] 


TJS  Laser 

EC  Laser 

</>ASE/(/>mcaS  *  22  °C 

0.17 

0.01 

</>ASE/</>mMsat32°C 

0.31 

0.02 

</>ASE/</)meas  at  38  °C 

0.49 

0.05 

</)ase/</>-«.  at  46  °C 

0.79 

0.38a 

the  vapor  when  the  main  lasing  mode  is  tuned  on  resonance. 
In  our  measurements  of  the  RIN,  we  were  interested  in  the 
ratio  of  the  intensity  fluctuations  to  the  intensity  of  the  main 
laser  mode,  which  was  the  mode  interacting  with  the  atoms 
in  the  vapor.  Thus  defining  (/)  Ase  as  an  estimate  of  the  la¬ 
ser’s  amplified  spontaneous  emission,  we  employed  {/) 
—  )meas~  (I)  ase  in  the  definition  of  the  RIN  noted  above. 
We  estimated  the  magnitude  of  the  ASE  by  increasing  the 
vapor  cell  temperature  to  56  °C,  and  then  measuring  the 
transmitted  intensity  of  each  laser.  Using  the  optical  depth 
measurements  from  the  lower  temperatures,  we  were  able  to 
infer  that  the  optical  depth  at  56  °C  was  approximately  1  cm, 
short  enough  so  that  no  intensity  from  either  laser’s  main 
mode  should  have  passed  through  the  vapor.  We  found  that 
the  TJS  laser’s  amplified  spontaneous  emission  was  about 
11%  of  the  total  laser’s  output  [16],  while  the  EC  laser’s 
ASE  was  only  about  1%  of  the  total  output  intensity.  Table  I 
shows  the  on-resonance  values  of  (/)ASE^)meas  f°r  eac^  of 
the  lasers  at  the  four  temperatures  studied. 


HI.  RESULTS 

Figure  3  shows  the  measured  RIN  of  the  two  lasers  as  a 
function  of  laser  detuning  for  four  different  temperatures.  At 
low  temperatures,  where  the  vapor  is  optically  thin,  the  RIN 
displays  a  double-peaked  dependence  on  laser  detuning. 
Heuristically,  if  L(v-v0)  is  defined  as  the  absorption  line 
shape  centered  at  v0 ,  then  for  PM-to-AM  conversion  we 
expect 


(8I/(I))**{8vL)ms 


dL 

dv 


(1) 


where  (Sv^^  corresponds  to  the  laser’s  root-mean-square 
frequency  fluctuations.  In  words,  the  laser’s  phase  fluctua¬ 
tions  cause  the  laser  frequency  to  vary  stochastically  across 
the  atom’s  absorption  profile,  and  as  a  consequence  the 
amount  of  light  transmitted  by  the  vapor  becomes  a  stochas¬ 
tic  quantity.  The  absolute  value  of  the  absorption  line 
shape’s  derivative  implies  a  double-peaked  curve  for  PM- 
to-AM  conversion  RIN  as  a  function  of  laser  detuning,  and 
given  the  excited-state  hyperfine  structure  associated  with 
this  transition,  the  double-peaked  structure  could  be  expected 
to  exhibit  some  asymmetry.  The  important  point  of  this  ob¬ 
servation  is  that  neither  shot  noise  nor  the  laser’s  intrinsic 
RIN  would  be  expected  to  exhibit  such  a  double-peaked 
structure  [17].  Consequently,  the  double-peaked  appearance 
of  the  RIN  data  shown  in  Fig.  3,  along  with  the  inverse 
relationship  between  the  magnitude  of  the  RIN  and  the  va¬ 
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FIG.  3.  RIN  (multiplied  by  106)  vs  laser  detuning.  Black  dia¬ 
monds  correspond  to  the  results  obtained  with  the  EC  laser,  while 
circles  correspond  to  the  results  obtained  with  the  TJS  diode  laser. 

por’s  optical  depth,  provide  strong  evidence  that  our  mea¬ 
surements  of  the  RIN  are  indeed  probing  the  PM-to-AM  con¬ 
version  process. 

One  of  the  more  intriguing  aspects  of  Fig.  3  concerns  the 
relative  magnitude  of  the  RIN  for  the  two  lasers.  Given  that 
the  line  widths  of  the  two  lasers  differed  by  about  102  to  103, 
we  had  originally  anticipated  that  the  RIN  of  the  EC  laser 
would  be  insignificant  on  the  scale  of  the  RIN  for  the  TJS 
laser.  Specifically,  since  for  a  singlemode  laser  with 

a  nearly  Lorentzian  line  shape  scales  like  \jA  vL  [18],  all 
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FIG.  4.  Ratio  of  the  TJS  laser’s  RIN  to  that  of  the  EC  laser: 
open  triangles  correspond  to  7=38  °C,  while  filled  circles  corre¬ 
spond  to  7=  32  °C.  The  dashed  line  corresponds  to  the  anticipated 
value  for  this  ratio  based  on  an  EC  laser  linewidth  of  100  kHz. 


other  things  being  equal,  Eq.  (1)  implies 

<*//</»tjs  /AvIJS\,/2 

(*//</»EC  \A^F/ 


(2) 


At  a  minimum,  we  therefore  expected  this  ratio  to  be  10, 
and  quite  likely  to  be  nearer  to  25  given  the  EC  laser’s  speci¬ 
fied  100- kHz  linewidth.  Moreover,  this  ratio  was  expected  to 
be  independent  of  laser  detuning  from  resonance.  As  illus¬ 
trated  in  Fig.  4,  however,  when  we  examined  the  experimen¬ 
tal  ratio  of  the  RINTJS  to  RINec,  we  found  that  on  average 
this  ratio  was  ~6  for  the  lasers  tuned  on  resonance,  while  off 
resonance  this  ratio  dropped  to  — 3.  Thus,  while  the  effi¬ 
ciency  of  PM-to-AM  conversion  is  certainly  a  function  of 
laser  linewidth,  the  experimental  data  indicates  that  the  heu¬ 
ristic  description  of  PM-to-AM  conversion  embodied  by  Eq. 
(1)  overly  simplifies  the  fundamental  physical  process. 


IV.  COMPUTATION 

In  the  regime  of  weak  fields,  such  that  the  optical  transi¬ 
tion  is  not  saturated,  PM-to-AM  conversion  in  optically  thick 
vapors  can  be  described  semiclassically  [13].  Basically,  for  a 
three-level  atom  like  that  shown  in  Fig.  2(b),  where  optical 
pumping  may  be  operative,  laser  phase  fluctuations  cause 
fluctuations  in  the  real  and  imaginary  components  of  the 
density-matrix  coherence  cr3] .  Since  the  coherence  is  related 
to  the  vapor’s  electric  susceptibility  *,  these  laser-induced 
atomic  coherence  fluctuations  essentially  result  in  fluctua¬ 
tions  of  the  vapor’s  absorption  cross  section. 

Defining  u  and  v  as  the  real  and  imaginary  components  of 
cr31 ,  respectively,  and  cr,,-  sis  the  population  in  |i),  it  is  pos¬ 
sible  to  write  the  density-matrix  components  in  terms  of  their 
average  values,  and  laser-induced  fluctuations,  Sx(t): 


CTU=:'1+  S„(t), 

(3a) 

U  =  M0+  Su(t), 

(3b) 

v  =  v0+Sv(t). 

(3c) 

As  discussed  in  Ref.  [13],  the  average  values  are  given 

47hfs(rflrfc+A2) 

[8rhfs(rflrfc+A2)+rflfi2]’ 

2flAyhfs 

Uq  —  '  ~~  » 

[8rhfS(rarfr+A2)+rflfi2] 


by 

(4a) 


(4b) 


2fir  a^hts 

[8yhfs(rarfc+A2)+ran2]’ 


where 


A 

ra=rco,+ y +  r- 


(4c)  • 


(4d) 


a  a2 

rb-rcol+  2  +  2a  +  r; 


(4e) 


while  the  laser  induced  fluctuations  are  described  by  the  fol¬ 
lowing  set  of  stochastic  differential  equations: 


.  n<o 


(5a) 


Su  =  -  rcol+  -  <5U  +  yuo“(A+  8<o)8v-  8g)Vo,  (5b) 
1 


rcoi+T+ 


Sl(t)2\ 
2  A  I 


8V  +  yvG  +  (A  +  8o))8u  +  8u>Uq 


n(t) 

2  8°' 


(5c) 


In  these  expressions,  Tco\  is  the  collisional  dephasing  rate  of 
the  optical  transitions;  A  is  the  Einstein  A  coefficient  for  the 
excited  state;  <o0  is  the  average  laser  frequency,  while  80) 
represents  the  laser’s  stochastic  frequency  fluctuation;  2y 
represents  the  linewidth  of  the  laser  (as  will  be  discussed 
subsequently);  w31);  y^s  is  the  hyperfine  relax¬ 

ation  rate  between  |2)  and  |l);  and  ft  is  the  optical  Rabi 
frequency.  Note  that  in  Eqs.  (4)  the  average  value  of  the  Rabi 
frequency,  ft,  is  employed,  while  in  Eqs.  (5)  the  Rabi  fre¬ 
quency  is  allowed  to  take  on  the  characteristics  of  a  stochas¬ 
tically  varying  parameter. 

The  frequency  and  intrinsic  amplitude  fluctuations  of  the 
laser  are  modeled  as  Omstein-Uhlenbeck  processes  [19]. 
Specifically,  if  ft(z  =  0,f)  is  defined  as  the  Rabi  frequency  at 
the  entrance  to  the  resonance  cell  at  time  f,  then 

ft(0,f)  =  ft(0)[l +*(/)],  (6a) 


(x(t)x(t+T))=  —  exp(  — a>i|r|),  (6b) 

(8u)(t)8a)(t+T))=yP  exp(-/3|r|),  (6c) 
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FIG.  5.  Computational  results  of  the  RIN  vs  the  laser  linewidth 
parameter  2y.  All  is  the  optical  dephasing  rate  due  to  spontaneous 
emission,  while  A  vD  is  the  Doppler  linewidth  (FWHM). 

where  fi(0)  is  the  average  Rabi  frequency  at  the  resonance 
cell  entrance;  /5T1  is  the  correlation  time  of  laser  frequency 
fluctuations;  2 y  is  essentially  the  linewidth  of  the  laser  for 
our  nearly  Lorentzian  line  shapes  (i.e.,  y<P),  and  \!<ox  is 
the  correlation  time  of  the  laser’s  intrinsic  amplitude  fluctua¬ 
tions.  Note  that  in  this  model,  the  laser’s  intrinsic  RIN  in  a 
1-Hz  bandwidth  is  4  \[y/(ox . 

The  above  density-matrix  equations  are  a  function  of 
propagation  distance  into  the  vapor  cell,  z ,  and,  as  discussed 
more  fully  in  Ref.  [13],  they  are  solved  iteratively  to  propa¬ 
gate  the  electric  field  amplitude  through  the  vapor.  Briefly, 
for  a  small  increment  in  propagation  distance,  Sz ,  the  laser 
electric  field  E  can  be  expanded  in  a  Taylor  series: 

E(z+Sz)  =  E(z)-4kNi£l3((vo+6v))v8z9  (7) 

where  N  is  the  number  density  of  atoms  in  the  vapor,  ju,l3  is 
the  1 1 )  —  1 3 )  electric  dipole  moment,  k  is  the  wave  vector  of 
the  laser  light,  and  (*••)„  implies  a  velocity  average  over  the 
Maxwell-Boltzmann  distribution.  At  various  computational 
time  steps,  values  for  the  stochastic  laser  electric-field  am¬ 
plitude  and  frequency  at  z  =  0  are  determined  by  the  meth¬ 
odology  of  Camparo  and  Lambropoulos  [20].  The  propaga¬ 
tion  of  this  wave  front  through  the  medium  is  then 
determined  with  the  aid  of  Eqs.  (4),  (5),  and  (7).  Equations 
(5)  are  solved  at  each  position  z  by  a  fifth-order  Runge- 
Kutta-Fehlberg  method  with  adaptive  (temporal)  step  size 
[21],  and  the  step  size  is  restricted  so  that  it  is  always  smaller 
than  0.1A-1  and  0.ir~oJ.  Ten  sets  of  density  matrix  equa¬ 
tions  are  solved  at  each  position  z,  corresponding  to  specific 
velocity  subgroups.  The  attenuation  coefficient  (i.e.,  SEI Sz) 
is  averaged  over  the  Doppler  distribution  by  use  of  Gaussian 
quadrature  [22],  and  the  field  is  then  propagated  through  the 
atomic  vapor.  The  intensity  of  the  laser  at  the  exit  of  the 
medium  is  averaged  over  a  time  interval  equal  to  the  smaller 
of  either  A-1  or  r”0},  and  is  then  subjected  to  statistical 
analysis  to  determine  SI  in  a  1  -Hz  bandwidth  at  low  Fourier 
frequency. 

The  results  from  this  one  dimensional  model  are  shown  in 
Fig.  5,  where  the  computed  RIN  is  shown  as  a  function  of 
the  laser  linewidth  parameter  2y.  Other  parameters  used  in 
these  calculations  are  collected  in  Table  II,  where  it  is  to  be 
noted  that  rdIL  for  a  100-kHz  linewidth  laser  is  0.45.  The 
error  bars  shown  in  Fig.  5  correspond  to  our  uncertainty  in 


TABLE  II.  Parameters  employed  in  the  numerical  simulation. 


Parameter 

Value 

Laser  input  intensity 

1  /xW/cm2 

Resonance  cell  temperature 

40  °C 

Resonance  cell  length 

3  cm 

Rubidium  vapor  pressure 

2.9X  10~6  torr 

Estimated  Td  for  2y=0.1  MHz 

1.35  cm 

Residual  N2  buffer  gas  pressure 

10" 6  torr 

^  ^Doppler 

512  MHz 

^  ^Natural 

5.67  MHz 

Laser  intrinsic  RIN 

io- 10 

Laser  ft  parameter 

3  GHz 

Sz  for  computation 

0.03  cm 

estimating  the  spectral  density  of  amplitude  fluctuations  for 
the  transmitted  beam,  as  discussed  in  Ref.  [13]. 

Perhaps  the  most  striking  feature  of  Fig.  5  is  that  it  pre¬ 
dicts  that  the  RIN  is  a  nonmonotonic  function  of  laser  line- 
width.  Though  the  RIN  increases  like  \py  for  linewidths  less 
than  the  atoms’  intrinsic  dephasing  rate,  in  this  case  A/2,  for 
larger  linewidth  lasers  the  RIN  is  actually  a  decreasing  func¬ 
tion  of  y .  In  particular,  comparing  a  100-kHz  linewidth  laser 
with  a  60-MHz  linewidth  laser,  a  yfy  dependence  would  pre¬ 
dict  that  the  narrow  linewidth  laser  should  have  a  RIN  ap¬ 
proximately  25  times  smaller  than  that  of  the  broad  linewidth 
laser.  However,  the  data  of  Fig.  5  actually  indicate  that  a 
100-kHz  laser  may  exhibit  nearly  the  same  RIN  as  a  60-MHz 
laser,  in  reasonable  agreement  with  the  experimental  data  of 
Fig.  4,  where  the  ratio  of  the  RIN  values  on  resonance  is 
about  3  for  7=38  °C  (i.e.,  r^/L  =  0.53). 

In  order  to  pursue  the  origin  of  this  nonmonotonic  depen¬ 
dence,  it  is  necessary  to  recognize  that  the  RIN  depends  on  a 
ratio  of  noise  to  transmitted  light  intensity.  For  the  data  of 
Fig.  5,  these  are  displayed  separately  in  Fig.  6.  Figure  6(a) 
shows  the  average  transmitted  light  intensity  as  a  function  of 
laser  linewidth.  As  might  be  expected,  for  laser  linewidths 
less  than  the  homogeneous  linewidth,  optical  transmission 
increases  with  y  (i.e.,  the  peak  absorption  cross  section  de¬ 
creases  as  the  laser  phase  noise  adds  to  the  atom’s  total 
dephasing  rate).  However,  once  the  laser  linewidth  increases 
beyond  the  homogeneous  width,  it  interacts  with  a  greater 
number  of  velocity  subgroups,  leading  to  increased  absorp¬ 
tion  by  the  vapor.  Eventually,  when  the  laser  linewidth  be¬ 
comes  greater  than  the  Doppler  width,  the  transmitted  inten¬ 
sity  again  increases.  Note,  however,  that  the  RIN  is  inversely 
proportional  to  the  transmitted  light  intensity,  so  that  taken 
by  itself  Fig.  6(a)  would  predict  that  the  RIN  is  a  minimum 
for  linewidths  near  All.  Figure  6(b)  shows  the  noise  level  of 
the  transmitted  light  intensity  as  a  function  of  laser  line- 
width.  The  data  points  correspond  to  the  noise  values  used  in 
the  computation  of  the  RIN  for  Fig.  5;  the  solid  line  will  be 
discussed  subsequently.  Clearly,  the  striking  feature  of  Fig. 
6(b)  is  the  “fall-off’  of  FM-to-AM  conversion  noise  for 
laser  linewidths  greater  than  the  optical  dephasing  rate,  A/2. 
Together,  Figs.  6(a)  and  6(b)  indicate  that  the  nonmonotonic 
nature  of  the  RIN  for  linewidths  near  A/2  results  primarily 
from  a  nonmonotonic  relationship  between  the  efficiency  of 
PM-to-AM  conversion  noise  and  laser  linewidth.  Addition- 
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FIG.  6.  (a)  Filled  circles  correspond  to  the  computed  average 
transmitted  intensity  as  a  function  of  the  laser  linewidth  parameter, 
(b)  Circles  correspond  to  the  computed  transmitted  intensity  noise 
as  a  function  of  the  laser  linewidth  parameter,  while  the  solid  line 
corresponds  to  the  analytical  result  as  discussed  in  the  text. 

ally,  the  drastic  reduction  of  the  RIN  for  linewidths  greater 
than  the  Doppler  width  results  primarily  from  the  increased 
optical  transmission  by  the  vapor. 


V.  ANALYSIS 

In  this  section  we  consider  the  origin  of  the  nonmonotonic 
relationship  between  the  efficiency  of  PM-to-AM  conversion 
noise  and  laser  linewidth.  To  this  end,  we  note  that  PM- 
to-AM  conversion  derives  principally  from  fluctuations  in 
the  atomic  coherence  [13],  and  we  anticipate  from  Eq.  (7) 
that  the  stochastic  variations  of  transmitted  intensity  will  be 
proportional  to  the  stochastic  variations  in  the  coherence. 
Specifically,  assuming  that  the  coherence  is  approximately 
constant  along  the  propagation  path  at  any  instant,  Eq.  (7) 
yields 

l(L,t)=(l(L))exp[-SkNv.n(Sv(t))vL],  (8a) 

where 

(!(L))=I0exV[-SkNfLl3(v0)vLl  (8b) 


For  the  situations  of  practical  interest  (i.e.,  spectroscopy  and 
atomic  clocks),  we  can  further  assume  that  the  coherence 
fluctuations  are  relatively  small,  so  that  the  exponential  of 
Eq.  (8a)  may  be  expanded.  It  is  then  clear  that  SI(L,t) 
~(Sv(t))v. 


To  proceed,  we  restrict  our  consideration  to  fairly  weak 
fields  such  that  fl^y*^  implying  that  S  =  <5^=0,  and 

ra  =  rb^r.  Thus  defining  £(t)  as  ( Su  +  i8v )  and  £0  as 
(cr13),  Eqs.  (5)  result  in 

£=  —  [r  — i(A+  £<*>)]£+ (y+j£a>)£0.  (9) 

Ignoring  terms  of  second  order  in  the  fluctuating  parameters 
(e.g.,  So)£),  Eq.  (9)  becomes  a  complex  Langevin  equation: 

£+[T  — fA]£=(y+i<?fc>)£0,  (10) 

where  the  atomic  dephasing  rate  T  is  now  seen  to  act  as  a 
viscous  drag  on  the  atom’s  ability  to  follow  the  rapidly  fluc¬ 
tuating  stochastic  “force”  Sa)£ Q.  Integrating  Eq.  (10),  we 
obtain 


£(f)  =  £oe~(r~'A)'  l\y+iS<o)e{T-ili),'dt'.  (11) 

JO 


In  the  case  of  a  nearly  Lorentzian  line  shape,  where  we 
can  assume  that  /?  is  much  greater  than  the  atom’s  dephasing 
rate  or  laser  detuning,  the  laser  frequency  fluctuations  will 
appear  to  be  <5  correlated,  so  that  (Sa)(tf )  Sa)(t"))  =  2yS(t' 
-t").  Then,  defining  the  autocovariance  of  atomic  coher¬ 
ence  fluctuations 


C(t,t+ r)-<f(/)£*(f  +  r))-(m)(C*(t  +  r)>,  (12) 


it  is  straightforward  to  show  from  Eq.  (11)  that 


lim  C(t,t+  t)  =  C(t)  = 

t-+  °° 


fl2y 


i6r[(r+r)i+Ai] 


c-r|T|e-«Ar 


(13) 


Each  velocity  subgroup  in  the  vapor,  however,  will  have  its 
own  autocovariance  function  owing  to  the  different  Doppler 
shifts  of  the  atoms.  Averaging  Eq.  (13)  over  the  Doppler 
distribution  of  shifts  in  the  case  of  the  laser  tuned  to  the 
center  of  the  Doppler  broadened  line  shape,  and  assuming 
that  the  Doppler  width  is  much  greater  than  the  homoge¬ 
neous  linewidth,  yields 


<C(r)>„* 


f^yV'n*  ln(2) 

8r(r+y)AyD 


^-(2r  +  y)|r|< 


(14) 


Taking  the  Fourier  transform  of  (C(r))„,  we  obtain  the 
spectral  density  of  coherence  fluctuations,  S^,  which  is  pro¬ 
portional  to  the  spectral  density  of  intensity  fluctuations. 
Thus,  at  some  Fourier  frequency  /  within  a  small  bandwidth 
A /,  the  measured  intensity  noise  is  proportional  to 
where 


ty«)“  J_“jC(r)),C-’“Vr. 
Evaluating  the  Fourier  transform  yields 


f l2yyjTT  ln(2) 

(2T  +  y) 

4f(r+y)Av0 

(2  r+y)J  +  w5 

(15) 


(16) 


so  that 
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( SI  1(1) ) = 8  JTirkNfi  13L  /Sf( 2irf  ).  (17) 

Note  that  Eq.  (16)  predicts  that,  for  Fourier  frequencies  much 
less  than  the  atom’s  dephasing  rate,  the  spectral  density  of 
intensity  noise  will  be  white. 

The  solid  line  in  Fig.  6(b)  corresponds  to  the  intensity 
noise  predicted  by  Eq.  (16)  at  low  Fourier  frequencies,  ap¬ 
propriately  scaled  to  fit  the  computational  data.  Note  that  in 
the  regime  of  small  linewidths  (i.e.,  y^O,  Eq.  (16)  predicts 
that  the  noise  is  proportional  to  Vy,  while  in  the  regime  of 
very  broad  bandwidth  fields  (i.e.,  y^T)  the  noise  is  propor¬ 
tional  to  l/\/y.  Equation  (13)  provides  the  basic  rationale  for 
this  nonmonotonic  relationship  between  the  RIN  and  laser 
linewidth.  Though  the  variance  of  frequency  fluctuations  for 
a  singlemode  laser  (i.e.,  phase  diffusion  field)  scales  like  y, 
the  variance  of  atomic  coherence  fluctuations,  which  drive 
the  RIN,  has  a  nonlinear  dependence  on  y.  In  particular,  for 
the  homogeneous  case  of  Eq.  (13)  the  variance  of  atomic 
coherence  fluctuations  is  maximized  when  y=T. 

VI.  DISCUSSION 

We  have  examined  the  influence  of  single-mode  laser 
linewidth  on  PM-to-AM  conversion  experimentally,  compu¬ 
tationally  and  analytically.  The  significant  result  from  these 
studies  is  that  the  efficiency  of  PM-to-AM  conversion  is  a 
nonmonotonic  function  of  laser  linewidth.  Basically,  PM- 
to-AM  conversion  is  driven  by  atomic  coherence  fluctua¬ 
tions,  and  the  variance  of  these  fluctuations  (in  the  homoge¬ 
neous  case)  has  an  extremum  when  the  laser  linewidth  equals 
the  atomic  transition’s  homogeneous  linewidth. 

The  results  have  an  immediate  bearing  on  recent  attempts 
to  improve  the  performance  of  gas-cell  atomic  clocks  by  us¬ 
ing  diode  lasers  for  optical  pumping.  In  the  traditional  gas¬ 
cell  clock,  a  population  imbalance  between  the  ground-state 
hyperfine  sublevels  of  Rb87  is  achieved  by  optical  pumping 
with  a  rf  discharge  lamp  [23].  In  the  absence  of  a  microwave 
field,  resonant  with  the  ground-state  hyperfine  splitting,  the 
transmitted  intensity  of  the  discharge  lamp  is  at  a  maximum. 


However,  in  the  presence  of  resonant  microwaves  the  hyper¬ 
fine  population  imbalance  is  destroyed,  and  the  transmitted 
intensity  is  minimized.  Hence  the  intensity  of  transmitted 
light  in  this  device  provides  the  operative  signal  that  is  used 
to  stabilize  a  microwave  field  to  an  atomic  resonance. 

The  atomic  clock’s  short-term  performance  is  often 
couched  in  terms  of  the  signal -to-noise  ratio,  and  roughly  a 
decade  ago  it  was  predicted  that  more  efficient  optical  pump¬ 
ing  with  a  diode  laser  would  yield  orders-of-magnitude  im¬ 
provement  in  clock  stability  [24].  Though  experiments  have 
repeatedly  verified  the  superiority  of  diode  lasers  over  rf- 
discharge  lamps  with  regard  to  the  magnitude  of  an  atomic 
clock’s  signal,  they  have  also  indicated  that  PM-to-AM  con¬ 
version  noise  can  be  orders  of  magnitude  larger  with  stan¬ 
dard  diode  lasers  than  with  rf-discharge  lamps  [5].  Prior  to 
the  present  work,  this  result  was  completely  counter¬ 
intuitive,  as  the  linewidth  of  a  standard  diode  laser  is  roughly 
50  MHz  [11],  while  that  of  a  rf-discharge  lamp  is  -2  GHz 
[25].  However,  in  light  of  the  present  results,  and  the  fact 
that  optical  dephasing  rates  in  gas  cell  atomic  clocks  are 
—  100  MHz  [26],  it  can  be  appreciated  that  the  standard  di¬ 
ode  laser  devices  operate  in  a  regime  where  y~T,  the  worst 
possible  situation  for  PM-to-AM  conversion  noise.  Alterna¬ 
tively,  PM-to-AM  noise  is  considerably  reduced  by  employ¬ 
ing  a  rf-discharge  lamp  whose  linewidth  is  much  greater  than 
the  optical  dephasing  rate  [27].  Clearly,  for  improved  gas¬ 
cell  atomic  clock  performance  one  should  use  a  diode  laser 
with  a  narrowed  linewidth,  gaining  an  advantage  in  atomic 
clock  signal  amplitude,  while  paying  a  minimal  price  in  PM- 
to-AM  conversion  noise.  We  note  that  this  is  just  the  ap¬ 
proach  taken  by  Mileti  et  al. ,  who  demonstrated  a  short-term 
performance  for  the  gas-cell  atomic  clock  that  rivals  that  of 
the  passive  hydrogen  maser  [6]. 
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